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The flavor dependence of the naive time reversal odd ("T-odd") parton distributions 
for u- and d-quarks are explored in the spectator model. The flavor dependence of hi; 
is of significance for the analysis of the azimuthal cos(2<^>) asymmetries in unpolarized 
SIDIS and DY-processes, as well as for the overall physical understanding of the dis- 
tribution of transversely polarized quarks in unpolarized nucleons. As a by-product 
of the formalism, we calculate the chiral-odd but "T-even" function h 1L which enables 
us to present a prediction for the single spin asymmetry -A^J^ ' for a longitudinally 
polarized target in SIDIS. 

Naive time reversal-odd ("T-odd") transverse momentum dependent (TMD) parton dis- 
tributions (PDFs) have gained considerable attention in recent years. Theoretically they 
can account for non-trivial transverse spin and momentum correlations such as single spin 
asymmetries (SSA) in hard scattering processes when transverse momentum scales are on 
the order of that of quarks in hadrons, namely Pt ~ k± -C \J Q 2 , A prominent example is 
the Sivers function /f^ [TJ which explains the observed SSA in semi-inclusive deep inelastic 
scattering (SIDIS) for a transversely polarized proton target by the HERMES collaboration 
[2] . It describes correlations of the intrinsic quark transverse momentum and the transverse 
nucleon spin. The corresponding SSA on a deuteron target measured by COMPASS [3j van- 
ishes, indicating a flavor dependence of the Sivers function. Another leading twist "T-odd" 
parton distribution, the chiral-odd Boer-Mulders function hi; [4] correlates the transverse 
spin of a quark with its transverse momentum within the nucleon. We focus on the flavor 
dependence of these "T-odd" functions where for example hi; is important for the analysis 
of the azimuthal cos(20) asymmetry in unpolarized SIDIS and Drell-Yan 0[6]. We also 
consider the flavor dependence of the "T-even" function hi~ L , which is of interest in the 
transverse momentum and quark spin correlations in a longitudinally polarized target [7]. 

Considerable understanding of TMDs and fragmentation functions (FF) have been gained 
from model calculations using the spectator framework [8| \9\ \E\ l6| fT(J] fTT]. 

In this formalism we start (cf. [32]) from the definition of the unintegrated color gauge 
invariant quark-quark correlator which contains the gauge link indicated by the Wilson line, 
W[a | b] , and work in Feynman gauge in which the transverse Wilson line vanishes [13] . In the 
diquark model the sum over the complete set of intermediate on-shell states in the definition 
of the correlator is represented by a single one-particle diquark state | dq; pd q , A), where pd q 
is the diquark momentum and A its polarization. The diquark is "built" from two valence 
quarks which can be scalar-spin or axial vector-spin 1. The unintegrated correlator is then 
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Figure 1: Different vertices for the axial- vector diquark and contribution of the gauge link in the 
one-gluon approximation. 



given by 



8{{P- P ) 2 -m 2 s )G(P° -p°) 



0; P,S)=2^ — 1J (2 ^3 v — (P, S\ h (0) W[0 | oo, 0, Or] | dq; P - p, A) x 

(dg;P-p,A|W[oo ) 0,5r|0]^ i (0)|P,5). (1) 

The essence of the diquark spectator model is to calculate the matrix elements in Eq. JTJ 
by the introduction of effective The nucleon-diquark-quark vertices T S (N) and Y% X (N) are 
represented in Fig. [U(a). For example, the matrix element for the axial vector diquark is 

(adq; P-p;X\ ipi{0) \P, S) = i -=—s^(P - p; A) - — — , (2) 

V o P m q < 1" 

where the polarization vector of the axial- vector diquark is e M , u(P,S) denotes the nucleon 
spinor and M and m q are nucleon and quark masses, respectively. The unpolarized TMD /i 
is obtained by inserting these expressions into Eq. ([I]) and projecting from the quark-quark 
correlator 



h{x,f T ) = \jdp~ (Tr [ 7 +$(p;P,5)] +Tr [ 7 +$(p;P, -5)]) 



(3) 

-=xP+ 



We have also calculated the distribution of transversely polarized quarks in a longitudinally 
polarized target, hj L by replacing T + with r + rT 5 and S by Sl, the spin 4-vector in 
longitudinal direction. Their analytic expressions in the scalar and axial vector diquark 
sectors are given in |14j . 

In the spectator framework the "T-odd" TMDs [8 J are generated by the gauge link (HUH [6]. 
The leading contribution, arising from the interference between tree- and box graph which 
contains an imaginary part necessary for "T-odds", is represented in Fig. Q] (c) in which the 
double line is an eikonal, and I is the loop momentum. For an axial-vector diquark we model 
the composite nature of the diquark through an anomalous magnetic moment k [15]. In the 
notation of Fig. Q] (b) the gluon-diquark axial diquark vertex is 

= ~*e dq [g^ ( Pl + p 2 r + (1 + k) {g^ (p 2 + q)^ + 9^ (Pi - qf 2 )] ■ (4) 

For k = — 2 the vertex T ax reduces to the standard jWW -vertex. We express the matrix 
elements including the gauge link, (dq; P — p\ W[oo, 0, Ot | 0] ipi(Q) \P, S) in the one gluon 
approximation [14] . Projecting the Boer-Mulders function, 

M f 

e%p> T hi(x,&) = — dp' (Tr [$ unpo i(p, S)ia* + l5 ] + Tr [* unpo i(p, -S)ia t+ l5 ]) , (5) 
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where = e + ^ the axial-vector diquark contribution is given by the expression, 

~' ax (^. dL\ — 

8(2tt) 3 f T + m 2 P+ J (2tt) 4 \ 3 



4dkt-{*,&) = - : ( q :i g ' 2 ^ / I U (a + P ) 2 ) ^ (p 2 ) x 



2? pi; (P - p - ( £ 4(P - p; A)£ M (P - p; A) 



[g^ w • (2P - 2p - /) + (1 + k) (v a (P-p + l)P + vP(P-p- 2l) a )} 



[l-v + iO] [I 2 - A 2 + iO] [{I + p) 2 -m 2 + iO] 

r / \ 

Tr 



x 



(f + M) (V - R g ^P) W - m q ) 7 + 7 ,; {Hi + m q ) (V + R g ^j 75] } + h.c. . (6) 

V(P —p—l) denotes the propagator of the axial- vector diquark. Since the numerator in Eq. 
© contains at most the loop momentum to the fourth power we can write it in the following 
manner, J2Li N al- a J ai -l az + N (0) . The (real) coefficients (tensors) N^... ai depend only 
on external momenta and can be computed in a straight-forward manner. The integration 
over the light cone components, l + and l~, are easily performed; however, calculating the 
Z + -integral results in an integral that is potentially ill-defined. This happens when g(p 2 ) 
is a holomorphic function in p 2 and at least one of the Minkowski indices is light-like in 
the minus direction, e.g. ot\ = — , a>2,...,ati <E {+,_!_} resulting in an integral of the form 
j dl + S(l + )Q(—l + ), implying that l + = and l~ = 00. This signals the existence of a light 
cone divergence in Ref. [IT]. 

One can handle the light cone divergences by introducing phenomenological form factors 
with additional poles [14] . 

t 2 ,_ {p 2 ~rn 2 q )f{p 2 ) 
9aAP )~ ^p_ A 2 + i0 ]»- I 'J 

For n > 3 there are enough powers of l + to eliminate this divergence. f(p 2 ) is a covariant 
Gaussian [14] which cuts off the pt integrations and A is an arbitrary mass scale fixed by 
fitting /1 to data. Similarly, the Sivers-function is projected from the trace of the quark- 
quark correlator U) (see e.g. [16]). 

2S^ P yt T {x 1 f T ) = ^ / 'dp- (Tr [ 7 +$(p; P, S T )] - Tr [ 7 +$(p; P, -St)] ) 

1 J p+=xP+ 

(8) 

It is well-known [5] that in the scalar diquark approximation the and /f^ coincide. 
By contrast the different Dirac structure for the chiral even f^p and chiral odd in the 
axial-vector diquark sector, Eq. © and ([§]) respectively, lead to different coefficients in 
the decomposition Na}... ai [14] . We fix most of the model parameters such as masses and 
normalizations by comparing the model result for the unpolarized /1 for u and d quarks 
to the low-scale (/i 2 = 0.34GeV) data parameterization of GRV [17] see Fig. O Note that 
PDFs for u and d quarks are given by linear combinations of PDFs for an axial vector and 
scalar diquark, u = f/ sc + \ j ax and d = f ax [HCE]. For "T-odd" PDFs we fix the sign 
and the strength of the final state interactions, the product of the charges of the diquark 
and quark, by comparing for u and d quarks in the diquark model with the existing 

data parameterizations (see Ref. [18]). The "one-half" moments (where q — u,d) of 
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Figure 2: The unpolarized up and down quark distributions functions (left) versus x compared to 
the low scale parameterization of GRV [17] . The half moment of the Sivers function (center-left) 
and Boer-Mulders function (center-right) and the unpolarized up and down quark distributions 
(« = 1.0). The half-moment (left) and first moments (right) of the Boer-Mulders and Sivers 
functions versus x fitted to extractions from data were presented in Ref. [18] 



the Sivers and Boer-Mulders functions 



4 1/2) W 



i-L(?,l/2) 
J IT 
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y 2M JlT 



(x,p T 



are displayed along with the unpolarized u and d quark distributions in Fig. O The "one- 
half" and first moments [14] of the up and down quark Sivers functions are negative and 
positive respectively while the up and down quark Boer-Mulders functions are both negative 
over the full range in Bjorken-x. We also note that the u-quark Sivers function and Boer- 
Mulders function are nearly equal, even with the inclusion of the axial vector spectator 
diquark. 

Having explored the flavor dependence of the we are now in a position to extend early 
phenomenological work on "T-odd" contributions to azimuthal asymmetries in SIDIS[6]. We 
consider the spin independent double "T-odd" cos 2(f) asymmetry for tt + and n~ production. 
We focus on the important contributions to the cross section for unpolarized SIDIS [16] 



da 



2-na 2 



dx dy dz dtfih dP 2 ^ xyQ 2 



l-y + ^y) f uu,t + (1 - y) cos(2^) F™J 



COS 2(f>h 



where the structure function 2< ^ h involves a convolution of the Boer-Mulders and Collins 



fragmentation function 



pcos 2<ji h 

r uu 



2h ■ kxh ■ p T — kr • Vt 

MM h 



(9) 



C is the convolution integral. Our input for the Collins functions is based on recent work 
in [19j where the Collins function was calculated in the spectator framework. It was assumed 



that H^ dls iav ^ w _jj^-yJ av ) m pj on sec tor, thereby satisfying the Schaefer-Teryaev 
sum rule [20] locally. We estimate the azimuthal asymmetry A c ^ 2< ^ (cf. Eq. (J9j) ) , where 
A cos2<t> ^ j cog 2^ da/ J da. In Fig.lwe display the A™j} 2<A (P T ) in the range of HERMES [2] 
and future JLAB kinematics [21] as well as x dependence in the range 0.5 < Pt < 1-5 GeV/c. 
Having calculated the chiral-odd but "T-even" parton distribution h^ L we use this result to- 
gether with the result of Ref. [U] for the Collins function to give a prediction for the sin(2</>) 
single spin asymmetry Aul for a longitudinally polarized target. A decomposition into struc- 
ture functions of the cross section of semi-inclusive DIS for a longitudinally polarized target 



±(fav) 
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Figure 3: The cos 20 asymmetry for and ^ as a function of Pt at JLAB 12GeV and HERMES 
(center) kinematics. Right: cos 2<p asymmetry for as a function of x at JLAB 12GeV and 
HERMES kinematics. 



reads (see e.g. [16]) 

dauL 
dx dy dz d(f>h dP% ± 



2™ 2 
xyQ 2 11 



"(1 - y) sin(2^) F^ (20) + (2 - y)VT^sin(0, 



(10) 



| is the projection of the spin vector on the direction of the virtual photon. In a partonic 

picture the structure function F^ 2 ^ is a leading twist object (while is sub- leading) 

and given by a convolution of the TMD hj~ L and the Collins function (cf. |16j ) 



b UL - L 



2h ■ kx h ■ Pt ~ &t • Pt 

MM h 



(11) 



We display the results for the single spin asymmetry A s ^ 2 ^ in Fig. [4] using the kinematics of 
the upcoming JLab 12 GeV upgrade. We note that the n~ asymmetry is large and positive 



due to the model assumption H l 



-L(dis — f av) 



H 



. This asymmetry has been measured 



at HERMES for longitudinally polarized protons [22J and deuterons [23J. The data show that 
for the proton target at HERMES 27.5 GeV kinematics both 7r + and tt~ asymmetries are 
consistent with down to a sensitivity of about 0.01. These asymmetries could be non-zero, 
but with magnitudes less than 0.01 or 0.02. These results are considerably smaller than 
our predictions for the JLab upgrade. For the deuteron target the results are consistent 
with for 7r + and ir~ . This SIDIS data for polarized deuterons could reflect the near 
cancellation of u- and d-quark h± L functions and/or the large unfavored Collins function 
contributions. There is also CLAS preliminary data [24] at 5.7 GeV that shows slightly 
negative asymmetries for tt + and ir~ and leads to the extraction of a negative xhf^ . This 
suggests that the unfavored Collins function (for d —t n + ) is not contributing much. Data 
from the upgrade should help resolve these phenomenological questions. 

We have performed calculations of transverse momentum dependent parton distributions, 
including the Boer-Mulders function h^, the Sivers function along with h± L in the 
framework of an axial- vector and a scalar diquark spectator model. The calculation of these 
functions in both sectors allowed us to explore their flavor dependence, i.e. to compute hj~, 
f^r an d hj~ L for a u-quark and a ci-quark. We used these results along with the Collins 



and A 



fragmentation function H ^ . to estimate the azimuthal asymmetries A 
in SIDIS. In summary, a refined diquark spectator model, including axial vector di-quarks 



(cos(20)) 
UU 



(sin(20)) 
UL 
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Figure 4: Left Panel: The half-moment of xhf^ 1 ^ 2 ^ versus x compared to the unpolarized up and 
down quark distribution functions. Right Panel: The sin 2<j> asymmetry for ir + and ir~ as a function 
of x at JLAB 12GeV kinematics. 

leads to both u- and d-quark "T-odd" TMDs and provides the ingredients for predicting a 

range of asymmetries for future experiments. 
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